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Abstract 

Recent experiments revealed that the dielectric dispersion spectrum of fission 
yeast cells in a suspension was mainly composed of two sub-dispersions. The 
low-frequency sub-dispersion depended on the cell length, while the high- 
frequency one was independent of it. The cell shape effect was simulated by 
an ellipsoidal cell model but the comparison between theory and experiment 
was far from being satisfactory. Prompted by the discrepancy, we proposed 
the use of spectral representation to analyze more realistic cell models. We 
adopted a shell-spheroidal model to analyze the effects of the cell membrane. 
It is found that the dielectric property of the cell membrane has only a minor 
effect on the dispersion magnitude ratio and the characteristic frequency ratio. 
We further included the effect of rotation of dipole induced by an external 
electric field, and solved the dipole-rotation spheroidal model in the spectral 
representation. Good agreement between theory and experiment has been 
obtained. 
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I. INTRODUCTION 



The polarization of biological cells has a wide scope of practical applications like manip- 
ulation, trapping or separation of biological cells jl],3 , and thus an accurate characterization 
of the polarization is needed. While the polarization of biological cells can be investigated 
by the method of dielectric spectroscopy [[| as well as by the measurement of field-induced 
cell movements [|]|)|, the former method has a much higher resolution M. For biological 
cells, the main contribution to the dielectric dispersion is the structural (Maxwell- Wagner) 
polarization effects (T|. Because of the practical applications, there is a strong need for 
intuitive models as well as simplified equations which describe the parameter dependence of 
the polarization. Thus, various cell models have been proposed for the analysis of the po- 
larization mechanisms. However, due to the complexity of existing theories, these methods 
have not yet found broader acceptance. 

In this work, we propose the use of the spectral representation || for analyzing the cell 
models. The spectral representation is a rigorous mathematical formalism of the effective 
dielectric constant of a two-phase composite material ||. It offers the advantage of the 
separation of material parameters (namely the dielectric constant and conductivity) from 
the cell structure information, thus simplifying the study. From the spectral representation, 
one can readily derive the dielectric dispersion spectrum, with the dispersion strength as 
well as the characteristic frequency being explicitly expressed in terms of the structure 
parameters and the materials parameters of the cell suspension (see section II. B below). 
The actual shape of the real and imaginary parts of the permittivity over the relaxation 
region can be uniquely determined by the Debye relaxation spectrum, parametrized by the 
characteristic frequencies and the dispersion strengths. So, we can study the impact of these 
parameters on the dispersion spectrum directly. 

The plan of the paper is organized as follows. In the next section, we will review the 
spectral representation theory || and show that the dielectric dispersion spectrum of a 
cell suspension can be expressed in terms of the spectral representation. In section III, we 
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will apply the spectral representation to the various cell models and present an alternative 
approach. We show that a better agreement with the experimental data can be achieved. In 
section IV, we discuss the effects of dipole rotation. We will show that the dipole rotation 
effect has a strong impact on the dispersion spectrum when the cells are sufficiently long. 
Discussion on further applications of our theory will be given. 



We regard a cell suspension as a composite system consisting of biological cells of complex 
dielectric constant Z\ dispersed in a host medium of e^. A uniform electric field E = E z is 
applied along the z-axis. We briefly review the spectral representation theory of the effective 
dielectric constant to establish notations. 



where s = €2/ '(£2 denotes the relevant material parameter and rj(r) is the characteristic 
function of the cell structure. The electric potential 0(r) can be solved formally || 



where Gq(t — r') = l/47r|r — r'| is the free space Green's function. By denoting an operator 



II. FORMALISM 



A. Spectral representation 



The spectral representation is initiated by solving the differential equation 




(1) 




(2) 




(3) 



and the corresponding inner product 




(4) 
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it is easy to show that T is a Hermitean operator. Let s n and \n) be the eigenvalue and 

eigenfunction of V such that T|n) = s n \n), where < s n < 1 is a real eigenvalue. The 
integral equation can be solved symbolically: 

\z)E . (5) 



From the solution, we obtain the electric field and hence compute the effective dielectric 
constant e e in the spectral representation. We further define the reduced effective dielectric 
function ||: 

^(s) = i - r = -w<^>- (6) 

e 2 sv hiQ 

By inserting the complete set 1 = J2n \ n )( n \> we fi n d 



F n is defined as the spectral function: 

1 



v (z\n)(n\z), (8) 



which satisfies a sum rule ||: 



4E(^>(^> = i<«k> = = P, (9) 

where V\ is the total volume of the suspending cells and p the volume fraction of the cells. 

B. Dielectric dispersion spectrum 

For cells of arbitrary shape, the eigenvalue problem of the T operator can only be solved 
numerically. However, analytic solutions can be obtained for isolated spherical and ellip- 
soidal cells. For dilute suspensions of prolate spheroidal cells, the cells can be regarded as 
noninteracting. The problem is simplified to the calculation of s n and \n) with a single cell, 
which can be solved exactly. Only two of the F n are nonzero, due to the orthogonality of 
\n) with \z). 



Thus, in subsequent studies, we restrict ourselves to two poles (n = 1,2). From Eq.(^), 
the effective dielectric constant is written in the spectral representation: 

^(i-E^v <«>) 

\ n =l 6 

After substituting ei = e\ + ai/j2Trf and e 2 = e 2 + a 2 /j2nf into Eq. (|T0|) , where e and a 
are the real and imaginary parts of the complex dielectric constant, we rewrite the effective 
dielectric constant after simple manipulations: 

where e# and are the high-frequency dielectric constant and the low-frequency con- 
ductivity respectively, while Ae n are the dispersion magnitudes, f£ are the characteristic 
frequencies of the nth sub-dispersion. We obtain the dispersion magnitudes Ae n and the 
characteristic frequencies respectively 0: 

si(s-t) 2 



Aei = F x e 2 
Ae 2 = -F 2 e 2 



S(S - 8i)(t ~ Si) 2 ' 

s 2 (s-t) 2 
s(s - s 2 )(t - s 2 ) 2 ' 
a 2 s(t - Si) 



2iie 2 t(s — si) 
cr 2 s(t - s 2 ) 



2ne 2 t(s - s 2 ) 

where s = e 2 /(e 2 — ei) and i = <r 2 /(cr 2 — o"i). To compare with experiment data ||, we 

express the dispersion magnitude ratio and characteristic frequency ratio as 

A £l Fi s x (s - s 2 )(t - s 2 ) 2 



Ae 2 F 2 s 2 (s-si)(t-si) 2 ' 

n _ (t-s 2 ){ s - Sl ) 
n (t- Sl )(s- S2 y 



(12) 
(13) 



III. APPLICATIONS TO VARIOUS CELL MODELS 

In a recent work 0], we adopted the spheroidal model (SM) to analyze the cell suspen- 
sions. Here we briefly review the analytic results of the model: 

1 2 

si = L z , s 2 = L xy , Fx = -p, F 2 = -p. 



where 



q 2 - 1 (q 2 - I) 3 / 2 
(1 - L z )/2 




ln(g +v /g 2 -l), 



L 



are the depolarization factors along the z-axis and a;- (or y-) axis of the prolate spheroid and 
q is the ratio of length L to diameter D. 

In the spheroidal model, we neglected the presence of a cell membrane. To study this 
effect, we put forward the shell-spheroidal model (SSM) here. In this case, the biological cells 
are modelled as shell-spheroidal ones with a spheroidal core of complex dielectric constant 
ex, covered with a confocal spheroidal shell of e s . For a small volume fraction p of shelled 
spheroidal cells embedded in a host medium of complex dielectric constant e 2 , the effective 
dielectric constant e e is given by the dilute-limit expression: 



where b z is the dipole factor for a single-shelled spheroidal cell along the z-axis ||: 



where y is the volume ratio of core to the whole shelled spheroid, while b xy indicates the 
dipole factor along the x- (or y-) axis, which can be obtained by replacing the subscript z 
with xy in the expression of b z . As a matter of fact, the cell suspension consisting of shell- 
spheroidal cells dispersed in a host medium is a three-phase system. Although the spectral 
representation was generally valid for two-phase composites, we have recently shown that 
it applies to composites of coated spheres as well as to coated spherical particles randomly 
embedded in a host medium [|T^]. Note the sum rule J2 F n = p is no longer valid. Similarly, 
one can show that the spectral representaton also applies to the present system consisting 
of spheroidal cells with shells of complex dielectric constant e s dispersed in a host. The 
effective dielectric constant is then given by 



e e = e 2 + pe 2 (b z + 2b xy ). 




(e s - i 2 )[i s + L z {£i - e s )\ + (?i - e s )y[e s + L z (e 2 - e s )] 



3( 



e s - e x )(e 2 - e s )yL z (l - L z ) + [e s + L z {e x - e s )][e 2 + L z (e s - e 2 )] ' 




(14) 
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with N.P. being the nonresonant part which vanishes in the limit of unshelled spheroidal 
inclusions, where 

L z [l + (x - l)y + L z (-1 + x + y - xy)\ 



S2 = 

Fi = 
F 2 = 



x L z {x 


i) 2 (y- 


i) + Li(x-m y -iy 




-1)3/-^ 


- L xy (-l + x + y - xy)] 


X L X y {X — 


-mv- 


-l)+Ll y (x-my-iy 






px 2 y 


3[x — L z (x - 


- my 


-i)+Li(x-m y -i)r 






2px 2 y 



3[x - L xy (x - \y{y - 1) + LUx - 1)% - l)] 2 ' 



where x = e s /e2- We omit the complicated expression for the nonresonant part here. In 
what follows, for the sake of convenience, we assume: (1) y is a constant for all coated 
spheroid; (2) x is a real number. 

In Fig.l, we plot the structure parameters F n and s n versus x for various y and for (a) 
q = 3.46, (b) 7.17 and (c) 10.24, respectively. In all case, p = 0.01. We find F n is strongly- 
dependent on y for x > 0.5, whereas it is not the case for s n . It may be concluded that 
the dielectric property of the cell membrane has a minor effect on the dispersion magnitude 
ratio, but plays no role in the characteristic frequency ratio. 

To investigate the validity of these models, we compare to experimental data, which was 
extracted by using a temperature sensitive cell division cycle mutant of fission yeast, cdc25- 
22 ||. Asami's theory |J results are also plotted for comparison. From Fig.2, it is evident 
that our model gives a better comparison with experimental data than Asami's theory. The 
reason for the improvement lies in the introduction of the conductivity contrast t by using 
of the spectral representation. As stated in Ref. 0, the large difference between our model 
and Asami's theory is due to a large <j\ ^> a 2 used in contrast to Asami's claim o\ w o 2 |§. 
We further find that SSM provides a better fit than SM for the dispersion magnitude ratio 
Aei/Ae2, while SSM yields the same results as those of SM for the frequency ratio f^/fx 
(both curves overlap in the right panel of Fig.2), indicating that the dielectric property of a 
cell membrane is indeed unimportant. 
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IV. EFFECTS OF DIPOLE ROTATION 



According to the numerical results, we find that the SSM provides a better fitting with 
previous experimental data than SM, but this improvement is actually too small. In other 
words, the dielectric properties of a cell membrane does not play an important role in 
dielectric dispersion spectrum. But, those numerical result will also show that both SM and 
SSM cannot obtain a good agreement with experimental data. In the presence of an electric 
field, cells of large length may rotate in favor of the applied field, thus we propose another 
model, namely the dipole-rotation spheroidal model (DRSM) to obtain a better fitting. 

When the cells are long enough, the rotation of dipole becomes very important with the 
external electric field under consideration, and the system is in general anisotropic. We have 
to take into account the effect of dipole rotation on F\ and F2, even for a weak electric 
field. Let us compute them from a thermodynamic consideration. We will show that they 
in general depends on q. 

Consider a spheroidal cell in an electric field Eq. Its long axis makes an angle 9 with the 
field. The dipole energy of the cell is 



E d (q,6) = -Re 



q(b z cos 2 6 + b xy sin 2 1 



16 

where b z and b xy are dipole factors along and perpendicular to the long axis: 

1 , 1 



(15) 



2 s(L z -§y xy 3(L xy -§y 

Eq.flI5|) can be understood by the energy approach. For simplicity, suppose the major axes 
of the cells all lie along the electric field, i.e., 8 = 0, then the induced dipole moments of the 
cells give a contribution to the effective dielectric constant. In the dilute limit, 

i e = e 2 + 3pi 2 b z , 

where p = V\/V is the volume fraction of the cells. For a fixed external field condition, the 
total electrostatic energy density of the suspension is given by E t = — Re(e e i?o/87r), which 
is equal to — Re(e 2 -E , Q/87r) + E d /V, and hence the desired results. 
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We showed that the conductivity contrast t attains a small negative value and thus 
the complex material parameter s can be approximated by t. Consequently, both b z and b xy 
have positive values. The probability is given by the Boltzmann factor 

p(q,6) = Ae~ E ^ kBT (16) 

where A is a normalization factor such that / p(q,9)dQ = 1, where Q is the solid angle, 
dQ = sin 9d9dip. We can calculate F\ and F 2 by the following integrals 

= pj p(q, 9) cos 2 OdSl, F 2 (q) = pj p(q, 9) sin 2 9dtt. (17) 

The Fi(q)/F 2 {q) ratio may be obtained by integrating with respect to 9 from to 7r/2 by 
symmetry. In the absence of an electric field, E^(q, 9) = and p(q, 9) equals to a uniform 
distribution. In which case, we obtain Fi = p/3 and F 2 = 2p/3, and hence F\jF 2 = 0.5. 
When the electric field is weak enough, the ratio is still constant and F\jF 2 = 0.5. Otherwise, 
the ratio will increase rapidly with q. For q = 1, b z = b xy and Fi/F 2 = 0.5 always. The 
above result implies that both F± and F 2 depend strongly on q when there is an electric 
field. For large q, b z 3> b xy , the spheroids tend to align with the applied field and hence 
F\jF 2 becomes very large. 

It is found that the mean cell length depends on the cultivation time, whereas the 
diameter is almost unchanged in an experiment ||, which will be applied to compare the 
different models. In the following numerical calculation, without loss of generality, we neglect 
the small imaginary part of e 2 , and define a new parameter £: 

e - €2D3E ° (is) 

e "T6^T (18) 
which characterizes the electric field strength. 

We can readily obtain the dispersion magnitude ratio Aei/Ae2 as well as the dispersion 
frequency ratio by substituting the results of Fi(q) / F 2 (q) into Eqs.(|12"D and (|T3"D, and 

setting s\ = L z and s 2 = L xy . In Fig. 3, F\jF 2 is plotted versus q. It is shown that F\jF 2 
depends strongly on the axial ratio q, especially for large q or strong magnitude of external 
electric field. 



To compare with experimental data in Fig. 2, we obtain good agreement in the DRSM 
with £ = 0.017, which corresponds to a weak field E Q « 0.1 V/m. The results show that 
dipole rotation indeed plays an important role in the dielectric dispersion - we cannot neglect 
the effect of the rotation of dipole induced by the applied electric field, especially when the 
average length of cell is large. In addition, good agreement exists only for large cytoplasmic 
conductivity, as attributed to a higher ion concentration in their cytoplasm to avoid the 
shrinkage of cells due to a loss of water across the cell membrane. 



V. DISCUSSION AND CONCLUSION 

Here we would like to make a few comments. At low frequencies, the cell membrane 
effectively insulates the interior of cell. In other words, a potential builds up entirely over 
the cell membrane, leaving the interior of cell rather inactive to the field [IJ. Thus e\ <C e 2 
and s = 1 + and it is reasonable to use s = 1.001 as fitting parameter. On the other hand, 
we assume that the host medium has low loss and o~ 2 ~ and at the same time a large 
cytoplasmic conductivity o\ <t 2 . Thus t = 0~ and it is reasonable to use t = —0.0001 to 
fit the data. 

The resulting equations [Eqs. fll^) and (|13|)1 are indeed simple equations arising from 



the spectral representation. These equations serve as a basis which describe the parameter 
dependence of the polarization and thereby enhances the applicability of various cell models 
for the analysis of the polarization mechanisms. In this connection, the shell-spheroid model 
may readily be extended to multi-shell cell model. However, we believe that the multi-shell 
nature of the cell may have a minor effect on the dispersion magnitude ratio as well as on 
the characteristic frequency ratio. 

In the presence of external electric fields, field-induced motions such as rotation of cells, 
dielectrophoretic motion or vibrational motion may have a significant impact on the di- 
electric dispersion spectrum. With the recent advent of experimental techniques such as 
automated video analysis JTT as well as light scattering methods 0, the cell movements can 
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be accurately monitored. For purely rotational motions, the distribution of surface charge 
on the cell surfaces may deviate significantly from the equilibrium distribution for cells at 
rest, leading to a change in the polarization relaxation and in the dielectric dispersion spec- 
trum. In this regard, our recent work on dynamic electrorheological effects |12[], in which 
the suspended particles can have rotational motions, may be applied to cell rotational mo- 
tions. In Ref . , we found that the particles' rotational motions do change the polarization 
relaxation substantially. 

In this work, we considered a monodisperse cell suspension, in which the cells are of the 
same shape (i.e., same length and diameter). While the diameter of the cells may remain 
constant during the cultivation process, the cells may possess a wide distribution of cell 
lengths 0. A modified theory, which takes the distribution of length into account, is urgently 
needed and our spectral representation theory will certainly help. In this connection, we 
may apply a strong dc electric field (in addition to the ac probe field) to help separating 
the long cells from the short ones. Our results indicate that even in a moderate field, the 
long cells can easily be aligned with the applied field, while the short ones remain essentially 
randomly oriented. In this way, an emphasis of dispersion spectrum of the long cells can be 
made possible. 

In summary, prompted by the discrepancy between recent theory and experiment on 
fission yeast cells, we have proposed the use of spectral representation to analyze more 
realistic cell models. We adopted a shell-spheroidal model to analyze the effects of the 
cell membrane. It is found that the presence of a membrane has only a minor effect on 
the dispersion ratio, but plays no role in the frequency ratio. We further included the 
effect of rotation of dipole induced by an external electric field. It has been found that the 
dipole-rotation effect plays an important role in the dispersion magnitude, but it does not 
change the characteristic frequency ratio. We obtained good agreement between theory and 
experiment when dipole-rotation effect is included. 
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FIGURES 

FIG. 1. For SSM, F\, F2, s\ and s 2 are plotted versus the dielectric constant ratio x for different 
thickness parameter y: (a) q = 3.46; (b) q = 7.17; (c) q = 10.24. 

FIG. 2. Ratios of the dispersion magnitudes and the characteristic frequencies are plotted 
versus q. Asami's theory: o x k, a 2 ; SM: t = -0.0014, s = 5.0; SSM: t = -0.0014, s = 5.0, x = 2, 
y = 0.8; DRSM: t = -0.0001, s = 1.001, £ = 0.017 (i.e., E is about O.lV/m). Note that the 
curves of SM and SSM overlap in the right panel, while they are quite close in the left panel. 

FIG. 3. For DRSM, the ratio F1/F2 is plotted versus q for different electric field strength 
parameter £. 
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